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Spin accumulation caused by confining potential 
in a two dimensional electron system 
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We investigate spin accumulation caused by spin-orbit interactions (SOIs) in a two dimensional 
electron system confined in one direction. We calculate spin density caused by three kinds of SOIs; 
arising from edge potential, Rashba and Dresselhaus mechanisms. All SOIs are shown to generate 
out-of-plane spin accumulations at the lateral edges of the system. We also find that the Rashba and 
the Dresselhaus mechanisms are competitive. Especially, when the strength of both interactions are 
equal, their effects cancel out each other. In that case, only the edge potential mechanism becomes 
relevant, and analytical expression is obtained for the spin density. The edge potential mechanism is 
shown to induce spin accumulation similar to and consistent with the experiment [Sih et at, Nature 
Phys 1, 31 (2005)]. We discuss the mechanism of the spin accumulation for each SOIs in some 
detail. 

PACS numbers: 72.25.Dc, 73. 23. Ad, 85.75.-d. 
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I. INTRODUCTION 

Spin-orbit interaction (SOI) in semiconductors attracts 
much attention, since it makes possible to control elec- 
tron spins without magnetic field, a goal of spintronics 
[l], 0. SOI is very small in a vacuum, since its strength is 
inversely proportional to the Dirac gap creating electron- 
positron pairs ^ 10^ eV. In contrast to a vacuum, SOI is 
largely enchanced for electrons in semiconductors, since 
the corresponding energy of creating electron-hole paris 
is given by the band gap A ~ 1 eV Q. This is why SOIs 
play important roles in spintronics. 

SOIs have been shown to generate spin currents per- 
pendicular to bias current, a phenomenon called spin Hall 
effect (SHE), which is expected to be a useful tool for cre- 
ating spin polarized current. SHE is explained as a result 
of spin dependent sccaterings through impurity poten- 
tial by D'yakonov and Perel' [3], and independently by 
Hirsch [5] and Zhang Q. Meanwhile, other scenarios are 
proposed for p— type bulk semiconductors described by 
Luttinger Hamiltonian and for n— type two dimen- 
sional semiconductors with the Rashba type SOI ^S]. In 
addition to the Rashba SOI , Dresselhaus SOI [loj is 
also shown to generate spin current in exactly o ppo site 
direction to the one which the Rashba SOI gives [11| . 

These transverse spin currents have been expected to 
generate opposite out-of-plane spin accumulations at the 
lateral edges of sample, analogously to Hall effect; lon- 
gitudinal current induces charges with opposite signs at 
the transverse edges. This theoretical prediction was con- 
firmed in the recent experiments [H, [H, [H, [lB| . Several 
theoretical studies suggest that the spin accumulation are 




FIG. 1: Schematic figure of the system. A two dimensional 
system is sandwiched between two contacts, /in and /xl are 
chemical potentials and we assume /ih >Ml. 
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SOI [20| are shown to cause opposite spin accumulations 
at lateral edges. Recently, Xing et al. [1^ have shown 
that the SOI coming from the edge potential can induce 
the spin accumulation as well. Here, a question arises, as 
to how these three SOIs determine spin density. As far 
as authors are aware, no calculation including all three 
types of SOIs has been published. 



In this paper we study the effects of three SOIs on 
spin accumulation in a two dimensional electron system 
(2DES). Spin density is calculated analytically for the 
SOI from edge potential, and numerically for the Rashba 
and the Dresselhaus SOIs. Finally we investigate the ef- 
fects of all SOIs. All SOIs are shown to generate opposite 
out-of-plane spin accumulations at the lateral edges. We 
also find that the Rashba and the Dresselhaus SOIs are 
competitive, as is the case for SHE. The SOI from the 
edge potential is shown to generate similar acuumulation 
to the experiment [l3|- For each SOIs we discuss precisely 
the mechanism of spin accumulation. 
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FIG. 2: (a) The edge potential in real systems: Vroai(a;) is zero 
in the center and gradually increases toward the edges, (b) 
The edge potential in our study: considering only the edge 
effect 
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II. MODEL HAMILTONIAN 

We consider a ballistic 2DES with two reflcctionless 
contacts, illustrated in Fig[TJ It means that electrons 
injected from one contact move freely without scatter- 
ings in the 2DES, and enter the other contact without 
reflections. In our setup, electrons are confined in x (z) 
direction by edge potential K;dge(a;) {V{z)). Although 
the edge potential profile in a real system is like T4eai(2;) 
in FiglD^a), considering the edge potential like Knodci(a;) 
in Fig[2l^b) is enough to understand the spin accumula- 
tion observed in the experiments. It is because the spin 
accumulation occurs only in the vicinity of the edges, and 
as far as the edge spin polarization is concerned, there is 
no essential difference between these two potentials. We 
assume Vmodci(a;) to be harmonic potential to make cal- 
culation easier. In this model, Hamiltonian for electrons 
compose of free electron part Hq and the SOIs part iJsoi- 
Hamiltonian for free electrons is given by 
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Kiodci = -m*uj'^x^ 



(1) 

(2) 



where ~p is momentum, to* is effective mass of electrons 
(in our work, we choose GaAs as the host of the electronic 
syetem, i.e. to* = 0.067TOe), w is a harmonic oscillator 
frequency. We assume that the confining potential V{z) 
is so strong that all electrons are in the lowest energy 
state in z direction. Hamiltonian for the SOIs consists of 
three parts. 
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where H^dge, and Hd are respectively the SOI from 
the edge potential, the Rashba SOI, and the Dresselhaus 
SOI. (Ti (i = X, y, z) is Pauli matrix, g* is Lande factor 
{g* — —0.44 in GaAs), h is Planck constant divided by 
2tt, a and f3 characterize the strength of the Rashba and 



FIG. 3: Spin density and electron density are plotted as a 
function of x for Fermi energy Ef = 50 meV, hui — 5 meV, 
M-H — /iL = 1 meV. Up-spin (down-spin) electrons accumulate 
in — (+) X direction. I = s/hfrnFuj is the characteristic length 
of the harmonic oscillator. 



the Dresselhaus SOIs, A is 1.42 eV in GaAs. 
we take z axis as the spin quantization axis. 
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III. CALCULATIONS AND RESULTS 

In this section, we calculate spin density in z direction 
(ov) which was detected in the experiments [H, [H, HI, 
[l5j . To investigate the effects of each SOIs precisely, 
we divide this section into three subsections. Firstly we 
study the effect of only the SOI from the edge potential, 
and secondly the Rashba and the Dresselhaus SOIs, and 
finally, all of them. All calculations are done assuming 
zero temperature. 



A. SOI from the edge potential 

In the case of the SOI from the edge potential, Hamil- 
tonian is given. 



^ipl+Pl) + lm*u^x^+^^—^^a,PyX. (7) 



Here, we can replace cr^ and Py respectively, by s {+1 
for up-spin, —1 for down-spin) and hky, since these two 
operators commute with H. Before diagonalization, we 
rewrite Hamiltonian to see how the edge potential works 
on electrons through SOI, 
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The second term represents harmonic potential depend- 
ing on spin and wavenumber. Up- and down-spin elec- 
trons with the same wavenumber ky feel different po- 
tential effectively. For example, for up-spin (down-spin) 
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electron with —ky, its potential minimum shifts to — (+) 
X direction. This Hamiltonian can be easily diagonalized 
and eigenstates and energies are given by 



pikyV 



(9* - 1)^^ 
2m*A 



sky)xs , (9) 



En 



Huj(n 



(10) 

2m* 2 A " 

where $,1(2;) and Xs are wavefunctions of the harmonic 
oscillator for band index n and of spin part. As we can 
see easily, the center of density profile shifts depending 
on the electron spin in z direction and the wavenumber in 
y direction. Up- and down-spin states are energetically 
degenerate. Using these wavefunctions, we calculate spin 
density (az) summing up the expectation value for occu- 
pied state (n, ky,s), 



(11) 



An assumption of reflectionless contacts greatly simpli- 
fies the treatment of applied electric field: in FiglU elec- 
trons with positive (negative) momentum in y direc- 
tion occupy the states with energy lower than /zl (mh), 
since electrons injected from one contact go through the 
other side contact without reflections. In our assumption 
^J'H > ^J■h, electrons with energy < E < have neg- 
ative wavenumber, and contribute to spin accumulation 
at zero temperature. By assuming fin — /^l ^ Mh, /^l, we 
obtain analytical expression of spin density in the linear 
response regime, 

-\'^n,klin),s=-l{x)\'^){^J'H - A^l) , (12) 

where ky{n) is Fermi wavenumber of band index n given 

by 



Ef - hw{n+ i) 

(1 +1 (5* _ 1)2(^)2) 



(13) 



In this expression, \^ 



represents 



n,kl(n),s=+(-) 

that up-spin (down-spin) electron accumulates in — (+) x 
direction, so contributions from two states (n, fc^(n), s — 
±1) result in opposite spin accumulations. Of course, the 
summation for band index n gives the similar result. The 
spin density and the electron number density are plotted 
in FiglSl From this figure, it is seen that the SOI from 
the edge potential causes opposite spin accumulations at 
the lateral edges. The spin density profile depends on the 
Lande factor g* ^ 1. In the case > 1, spin density pro- 
file is opposite to the result in FigO We also plot the spin 
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FIG. 4: Spin density is plotted as a function of x and the edge 
potential strength hu for Ef = 50 meV, fiH — /iL = 1 meV. 
The width of the spin density and the numbers of oscillations 
depend on the trap strength. 
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FIG. 5: Spin density is plotted as a function of x and a for 
Ef = 50 meV, huj = 5 meV, hh ~ fiL ~ 1 meV, and (3 = 
5 X 10~^^ eVm. The Rashba and the Dresselhaus SOIs are 
competitive. Dash line indicates a = /3. 



density, as a function of the strength of the harmonic po- 
tential with keeping Fermi energy constant in Fig[4l The 
width of the spin density and the number of its oscilla- 
tions are characterized by the wavefunctions of electrons 
in the highest band. In the strong trapping limit, only 
electrons in the lowest subband generate the spin accu- 
mulation with the narrow width and a small number of 
oscillations. As the potential become weaker, the more 
subbands are involved and the spin density oscillations 
spreads over the wide region. The step like change oc- 
curs when the number of subbands crossing the Fermi 
level changes. Spin polarization, P = (az)/{n), at the 



peak of the spin accumulation is about 10 ^ 
FigOand FigH 
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B. Rashba SOI and Dresselhaus SOI 

We treat the Rashba and the Dresselhaus SOIs on an 
equal footing since they are theoretically shown to gen- 
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erate spin currents in opposite directions with the same 
magnitude, and expected to give similar effects on spin 
density. In our model, Hamiltonian is given by, 

+ -^{PxCTy - Pya^) + -{Pa;ax - PyfJy) , (14) 

Similar to the previous section, we can put py — hky. 
Wavefunctions is expanded in eigenfunctions of the har- 
monic oscillator, 

^N,k,M^,y)= E «r;*^*n,fe.,. , (15) 

where iV is a new band index, ± indicate energy bands 
split caused by the SOIs. Coefficients deter- 
mined numerically. We expect non-zero spin density in 
applied electric field as in the previous section. A cal- 
culation for spin density (cr^) is the same as Eq(ll). To 
see the effects of two SOIs, the spin density is plotted in 
FigO as a function of the Rashba strength with keep- 
ing the Dresselhaus strength constant, since the Rashba 
strength can be tuned by gate voltage [l^l- From FiglSl 
two SOIs are shown to generate competitive result, that 
for a > (3 {a < f3), the spin density at the edges is 
up-spin-like (down-spin-like) in +x direction, and down- 
spin-like (up-spin-like) in —x direction. In the case 
a = f], the effects of two SOIs cancel out each other 
and the spin density is exactly zero in the entire region. 
Spin polarization at the peak is about 10~^. This accu- 
mulation is obviously caused by two SOIs, but how the 
SOIs work physically is unclear. Such a difhculty comes 
from the fact that electron spins are not conserved. To 
see how the SOIs change electron spin states, we derive 
dynamical force acting on electrons in x direction, given 
by = m*(Px/dH = m* {ih)-^[[x, H], H] [H], 

p dVEdgc{x) 2m* 2 fl2M, 

= ^ ^ ~ ^ '^'^^ ■ ^ ^ 

The second term shows that the Rashba and the Dressel- 
haus SOIs cancel out in the case |a| = |/9|. In quantum 
mechanics, the force is hard to deal with, so we calculate 
back an effective potential from F^, just as in classical 
mechanics, giving a spin dependent effective potential. 



V,s{x) = - 




= VEdgo(a;) ^{a'^ - P'^)kyazX . (17) 

From this expression, it is clear that the second term 
determins the profile of spin density. 

We now consider energy bands to discuss the effect 
of the second term, as illustrated in Fig[6l A state 
A with energy £'n,fc„,s=±ij an eigenvalue of Hamilto- 
nian Hq, is degenerate for up- and down-spin. Spin 
density, contributed from the states {n,ky,s = 1) and 
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FIG. 6: Schematic figure of energy separation. The degener- 
ate state A with energy En^ky,±i, eigenenrgy of Ho, is separate 
into the stete B with higher energy E^ ,^^ j^ and the state C 

with lower energy E„ ,^^ _. D and G are the states with Fermi 
energy. 
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FIG. 7: Electron densities and spin densities are plotted as 
a function of x for state (10, fc^ , and (10, fc^ , — ) in highest 
two bands crossing Fermi energy. Parameters are E-p = 50 
meV, fttj = 5 meV, a = 7 x 10"^^ eVm, /3 = 5 x 10"^^ 
eVm. Solid line represents that the profiles of electron number 
densities are approximately same. Dashed line and dotted 
line show the doformed spin densities for two state caused by 
SOIs. 



(n,ky,s = —1), is exactly zero, since spin densities for 
up- and down-spin have the same profile with opposite 
signs, and cancel out each other. The SOIs lift the de- 
generacy of two states into a state B with higher enegy 
Ej^i j,^ _|_ and a state C with lower energy Ejq ky -■ 
the case ky<0 and |q;|>|/3|, for the state B (C) , the spin 
density, following the second term, behaves up-spin-like 
in a; > (x < 0) and down-spin- like in 2: < (x > 0), 
while the electron density is approximately unchanged. 
In Figl?! we plot the spin densities and the electron num- 
ber densities for two states, say states D and G. Here the 
states D and G are taken to be the two highest band 
states at the Fermi level. In our parameter, they corre- 
spond to the states labeled by (10, ky, +) and (10, ky, —). 
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Profiles of electron densities for two states are nealy the 
same (solid line), but profiles of spin densities are differ- 
ent as explained above (dashed and dotted line). In our 
system, the density of states is inversely proportional to 
the square root of energy measured from the bottom of 
subband. With the same band index N and energy E, 
electrons in + band contribute to spin density more than 
electrons in — band do, then the result for |q;|>|/3| in 
Figl5]is similar to the spin density profile for electrons in 
+ bands. In this way, we can explain the induced the spin 
accumulation with the discussion from the deformation 
of the spin density caused by the SOIs and the density 
of states. 



C. All SOIs 

In this subsection, we consider all SOIs. Hamiltonian 
is exactly diagonalized numerically by expanding wave- 
functions in eigenstates of the harmonic oscillator. Cal- 
culated spin density (tTz) is plotted in FiglSl representing 
the effects of three SOIs. This figure looks similar to the 
previous result (FigE]), but, even in the case \a\ = \(3\, 
the spin density is not zero. It is because the SOI from 
the edge potential makes spin density at the edges up- 
spin-like (down-spin-like) in — (-I-) x direction. In that 
case, the spin density profile is the exactly same result of 
only the SOI from the edge potential, as long as |a| — \/3\. 
We believe that a result obtained in the realistic poten- 
tial like Fig[2fa) is essentially same with our results, since 
only the physical mecanism at the edges is important. 

We obtain similar spin accumulation to the experiment 
(a = 1.8 X 10^^^ cVm, (3 — eVm), performed in the two 
dimensional electron system [14]. In that experiment, the 
observed spin density is up-spin-like (down-spin-like) at 
the edges in — (-I-) x direction. This shows the impor- 
tance of the SOI from the edge potential, since the theory 
taking account of the Rashba SOI alone can not explain 
the experimental result 14]. 

One of the most characteristic results is that the signs 
of spin accumulation from the edge potential is opposite 
depending on the Lande factor g* ^ 1. It is because the 
sign of Hedge (4) depends on the Lande factor. Actually, 
in the experiment performed in ZnSe {g* = 1.1) [isj . the 
sign of the spin accumulation is opposite to the exper- 
iment peformed in GaAs {g* = -0.44) [H, [ll|- These 
experimental results are consistent with our analysis, and 
demonstrate the important role of the SOI from the edge 
potential. 
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FIG. 8: Spin density caused by three SOIs is plotted as a 
function of x and a for Ef = 50 meV, hui — 5 meV, jJ.H—fJ'L = 
1 meV, and /9 = 5 x lO^^'^ eVm. Spin density is similar to 
the previous result in Fig[5] but for \a\ = \(3\ the spin density 
is not zero because of the SOI from the edge potential. 

that the Rashba and the Dresselhaus SOIs produce com- 
petitive results. Especially when the strength of both 
SOIs are equal, their effects cancel each other. In that 
case, the spin accumulation caused by the SOI from the 
edge potential remains, and the analytical expression was 
obtained for the spin density. In addition, we discussed 
precisely the physical mechanism of forming opposite 
spin accumulations for each SOIs, considering effective 
potentials acting on electrons. Finally, we discussed our 
results, compared with the experiments, and concluded 
that the SOI from the edge potential induces the observed 
spin accumulation. From our result, we predict that the 
spin accumulation for 17* < 1 {g* > 1) can be diminished 
(enhanced) as the Rashba strength increase, which can 
be observed by tuning the gate voltage. This electrical 
control of electron spins would become a useful tool in 
spintronics. 
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IV. SUMMARY 



We have investigated the spin accumulation caused by 
three SOIs in the 2DES. The spin density (cr^) was cal- 
culated including all SOIs. All three SOIs were shown to 
induce opposite spin accumulations at the lateral edges, 
which was observed in the experiments . We also showed 



6 



[1] G.A Prinz, Science 282, 1660 (1998). 

[2] S.A Wolf, D.D.Awschalom, R.A.Buhrman, 

J.M.Daughton, S.von Molna'r, M.L.Roukes, 

A.Y.Chtchelkanova, and D.M.Treger, Science 294, 

1488 (2001). 
[3] E. I. Rashba , Nature Phys. 2, 149 (2001). 
[4] M. I. D'yakonov and V. I. Perel', JETP Lett. 13, 467 

(1971). 

[5] J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 

[6] S. Zhang, Phys. Rev. Lett. 85, 393 (2000). 

[7] S. Murakami and N. Nagaosa, Science. 301, 1348 (2003). 

[8] J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jung- 

wirth, and A. H. McDonald, Phys. Rev. Lett. 92, 126603 

(2004). 

[9] Y. A. Bychokov and E. I. Rashba, J. Phys. C 17, 6039 
(1984). ' 

[10] G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

[11] N. A. Sinitsyn, E. M. Hankiewicz, W. Teizer, and J. 

Sinova, Phys. Rev B. 70, 081312 (2004). 
[12] Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 

Awschalom, Science. 306, 1910 (2004). 
[13] J. Wnderlich, B. Kaestner, J. Sinova, and T. Jungwirth, 

Phys. Rev. Lett. 94, 047204 (2005). 
[14] V. Sih, R. C. Myers, Y. K. Kato, W. H. Lau, A. C. Gos- 



sard, and D. D. Awschalom, Nature. Phys. 1, 31 (2005). 
[15] N. P. Stern, S. Ghosh, G. Xiang, M. Zhu, N. Samarth, 
and D. D. Awschalom, Phys. Rev. Lett. 97, 126603 
(2006). 

[16] M. Governale and U. Ziilicke, Phys. Rev B 66, 073311 
(2002). 

[17] B. K. Nikolic, S. Souma, L. P. Zarbo, and J. Sinova, Phys. 

Rev. Lett. 95, 046601 (2005). 
[18] Q. Wang, L. Sheng, and C. S. Ting |cond-mat /0505576J 
[19] J. Wang, K. S. Chan, and D. Y. Xing, Phys. Rev B 73, 

033316 (2006). 

[20] A. G. Mal'shukov, L. Y. Wang, C. S. Chu, and K. A. 

Chao, Phys. Rev. Lett. 95, 146601 (2005). 
[21] G. Usaj and C. A. Balseiro, EuroPhys. Lett. 72, 631 

(2005). 

[22] J. Yao and Z. Q. Yang, Phys. Rev B 73, 033314 (2005). 
[23] Y. Xing, Q. Sun, L. Tang, and J. Hu, Phys. Rev B 74, 
155313 (2006). 

[24] J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 

Rev. Lett. 78, 1335 (1997). 
[25] B. K. Nikohc, L. P. Zarbo, and S. Welack, Phys. Rev B 

72, 075335 (2005). 



